Gravitational waves from relativistic neutron star mergers with microphysical 

equations of state 
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The gravitational wave (GW) emission from a set of relativistic neutron star (NS) merger sim- 
ulations is analysed and characteristic signal features are identified. The distinct peak in the GW 
energy spectrum that is associated with the formation of a hypermassive merger remnant has a 
frequency that depends strongly on the properties of the nuclear equation of state (EoS) and on 
the total mass of the binary system, whereas the mass ratio and the NS spins have a weak influ- 
ence. If the total mass can be determined from the inspiral chirp signal, the peak frequency of the 
post-merger signal is a sensitive indicator of the EoS. 

PACS numbers: 95.85.Sz, 97.60.Jd, 95.30.Lz 



Among the strongest known sources of gravitational 
wave (GW) emission are the merging events of double 
neutron star (DNS) binaries. Recent population systhe- 
sis studies (e.g. [jy]) and the discovery of the DNS J0737- 
3039 suggest a possible detection rate of GW radia- 
tion from DNS mergers of one in ^^30 years for LIGO I 
and one every two days for advanced LIGO. To detect 
such GW signals and to filter them out of the detec- 
tor output, theoretical waveform templates are needed. 
While the inspiral phase prior to the actual merger can 
be described very accurately within the post-Newtonian 
(PN) framework (e.g. Q), hydrodynamical simulations 
are needed to model the dynamical merging phase. In 
addition, different aspects of physics enter the problem 
at this stage. Besides general relativity (GR), nuclear 
and particle physics play a role in the description of the 
hot and dense NS fluid via an equation of state (EoS) 
and in the treatment of energy losses (e.g., by neutrinos) 
after the merging. The GW signal of the late inspiral and 
merging phases is therefore expected to contain informa- 
tion not only on the binary parameters such as masses 
and spins but also on the nuclear EoS. 

Efforts to investigate NS mergers have concentrated ei- 
ther on the relativistic aspects while simplifying the mi- 
crophysics (e.g. Q and refs. therein), or have employed a 
microphysical EoS together with an approximative neu- 
trino treatment while describing gravity in a Newtonian 
framework (e.g., The conformal flatness approach, 

a middle ground between PN and full GR, combined with 
a nuclear physics-based nonzero-temperature (T ^ 0) 
EoS has recently been chosen by Oechslin et al. 0] ■ 

The generic GW signal from a NS merger can be split 
into a chirp-like part emitted by the inspiralling binary, 
the burst amplitude from the final plunge when the two 
stars collide (when time is set to t = in Fig. [1]) , and 
a quasi-periodic post-merger signal caused either by the 
rotation and internal oscillation of a newly formed, non- 
axisymmetric hypermassive NS (HMNS) as merger rem- 
nant, or by the quasinormal ringing of a newly born black 
hole (BH) in case of a prompt gravitational collapse of 



the remnant after the final plunge. A first maximum of 
the compactness of the relic HMNS is associated with 
a minimum of the amplitude h = (|/i+P + |/ixP)^^^ 
at about 0.5 ms after the merging, followed by the on- 
set of the characteristically different, quasi-periodic post- 
merger emission. For some of our computed models, the 
quantity h is plotted in Fig. [1] It contains the com- 
bined information from both polarisations and hx oi 
the wave amplitude and therefore represents the enve- 
lope of the high-frequency wave pattern. Its post-merger 
modulation is caused by the oscillation of the nonaxisym- 
metric remnant. Since the pre- and post-merger signals 
are emitted in different frequency bands, they can be 
clearly identified in the corresponding luminosity spec- 
trum. The inspiral signal leads to a broadband contribu- 
tion below '^l kHz and depends mainly on the NS masses 
and their spins, while an EoS dependence is only present 
in the very last stage before merging fs*]. On the other 
hand, the post-merger signal is dominated by a quasi- 
periodic wave pattern with a frequency of about 2-4 kHz 
or about 6-7 kHz, depending on whether a HMNS forms 
or a prompt collapse to a BH happens 0. The associ- 
ated peak in the luminosity spectrum can become very 
pronounced in cases where the remnant keeps radiating 
GWs for several tens of ms as suggested by recent merger 
simulations 0, 0| • The bare presence of a contribution in 
the frequency range of about 2-4 kHz indicates the for- 
mation of a HMNS and a nuclear EoS that is sufficiently 
stiff to prevent prompt BH formation. 

In the present letter, we concentrate on the HMNS 
formation case and assess the question, to which ex- 
tent the nuclear EoS and the binary parameters can be 
constrained when such a post-merger peak is detected. 
Based on a set of simulated binary merger models [tI, 
we identify characteristic features of the simulated GW 
signals and link them to the merger properties. The sim- 
ulations were carried out with our relativistic smoothed 
particle hydrodynamics (SPH) code l3, U | , which solves 
the relativistic hydrodynamics equations together with 
the Einstein field equation in the conformally flat ap- 
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FIG. 1: The GW amplitude h for difFerent EoSs and spins as 
radiated perpendicularly to the orbital plane of the merging 
binary and measurable at a distance of 20Mpc. The mini- 
mum at about 0.5 ms is considered to mark the onset of the 
quasi-periodic wave train caused by oscillations of the rapidly 
spinning, non-axisymmetric merger remnant. 



proximation [CFC;[l3, The simulations were started 
from a stable equilibrium configuration slightly outside 
the innermost stable circular orbit and the correspond- 
ing initial data were generated by relaxing the fluid to 
a velocity field that includes the orbital motion and the 
proper spins of the NSs. Two T 7^ EoSs, the Shen-EoS 
[3], and the Lattimer-Swesty-EoS an ideal-gas EoS 
with parameters chosen to mimic the Shen-EoS, and the 
APR-EoS ^ were used. The APR-EoS was extended 
by an ideal-gas-like thermal pressure contribution that is 
proportional to the internal energy increase due to shock 
heating and viscous heating [Ij. The size of this contri- 
bution is determined by an adiabatic index Pth for which 
we chose two different values (Fth — 1.5, 2) in order to 
investigate its influence on the merger outcome. Finally 
we calculated two models with the Shen-EoS, restricting 
the latter to T = in order to investigate the influence 
of temperature-dependent pressure terms (see [7] ) . 

The Shen-EoS is relatively stiff and for NSs with 
typical masses around 1.4 M© leads to radii as big as 
^ns ^14 km. Its maximum mass of non- rotating NSs is 
~2.25Mq. In contrast, the LS-EoS is much softer and 
yields a radius near 12 km for a 1.4 M© NS. The APR 
EoS is still softer below and around nuclear density, but 
becomes very stiff at higher densities 3 x 10"'^''gcm~'^). 
Therefore it makes NSs even more compact (i?ns ~ 
11km), although it allows for a rather large maximum 
NS mass of ^^2.2^/© compared to about 1.8 Af© for the 
LS-EoS (see Fig. 2 in [7]). Besides the EoS, we have 
also varied the NS masses, the mass ratio, and the NS 
spins in our calculated set of models (see Table H]). 

The GW waveform hij is extracted by making use 
of the quadrupole formula and is given by hij = 
{2/ D)cP'Qij / dt^ ^ where Qij is the Newtonian mass 



TABLE I: Characteristic quantities of our computed models. 
Models with names starting with 'S' use the T 7^ Shen- 
EoS, 'C' models the restriction of this EoS to T = 0, the 'LS' 
model uses the LS-EoS, the 'P' model the ideal-gas EoS, and 
the 'A15' and 'A2' models the APR-EoS extended by ideal 
gases with different values of Fth . All models were computed 
with irrotating initial conditions except the last four cases 
where the ending 'co' ('ct') of the model names indicates ini- 
tially corotating (counterrotating) spin states of the NSs (for 
the spin frequencies, see Q)- Afi and M2 are the individ- 
ual gravitational masses in isolation, and q — M1/M2 is the 
mass ratio, /max, /peak, A_Bin, and AEpm are defined in the 
text. 'SNR' means the estimated signal-to-noise ratios in ad- 
vanced LIGO (left) and DUAL (right) for the GW emission 
after merging and a source distance of 20 Mpc. 
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quadrupole and D is the distance from the source (the 
indices ij denote the spatial directions). Compared to a 
more detailed extraction technique in the wave zone us- 
ing the gauge- invariant Moncrief variables as done, e.g., 
in [4] , this approximation is able to describe the GW sig- 
nal only qualitatively. The wave phase and thus the fre- 
quency information can be well reproduced but the am- 
plitudes are underestimated by about 30% in the inspiral 
regime and by about 40% in the post-merger regime [13, 
Fig. 12, panel a]. Based on the thus obtained waveform, 
the GW luminosity spectrum can be determined accord- 
ing to [li by d£;Gw/rf/ = lATrD^P{\h+\^ + \hx\^), where 
/i+ and hx denote the fourier transforms of the wave- 
forms of '+' and 'x' polarisation, respectively. The an- 
gle brackets indicate averaging over all possible source 
detection angles. The energy emitted in GWs is then 
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FIG. 2: GW luminosity spectra for different total masses, 
mass ratios, spin configurations, and EoSs (from top to bot- 
tom) . Note that the spectrum below 1 kHz is not represented 
correctly because we started the simulations when the GW 
frequency was ~ 1 klfz. 



given by Ai?GW = / df dEQy^/df. Because of the un- 
derestimation of the GW amplitude, the GW luminosity 
spectra, which depend quadratically on the amplitude, 
are systematically too low by ~70%. 

Bearing this in mind, we consider in the following 
quantities that are not directly affected by this shortcom- 
ing and independent of the source orientation, namely 
(see Table HI: (i) /max as the frequency of the GW signal 
when the amplitude becomes maximal at about the time 
of the final plunge. It is determined by fitting a function 
of the form A{t) cos{Lo{t)t + (f)) to the waveform; (ii) /peak 
as the frequency of the post-merger peak in the luminos- 
ity spectrum; (iii) the ratio AE'in/AE'pm of AE^in as the 
energy emitted over a time interval of 3 ms before merg- 
ing, and A_Epni as the energy radiated over a time interval 
of 5 ms after merging. The energies are determined as de- 
scribed above from the waveforms produced in the corre- 
sponding time intervals. The values thus obtained agree 
with time integrals of the quadrupole-formula-based ex- 
pression dE/dt = l/5(QijQij) to within ~20%. 

In cases where a HMNS forms, /peak turns out to de- 
pend sensitively on the EoS (Fig. [21 bottom) and to 
a lesser extent on the total mass of the binary system 
(Fig. [21 top) . The NS spins and the mass ratio have very 
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FIG. 3; Ratio Ai5pm/Ai5in vs. /peak (top) and /m ax VS . /p eak 
(bottom) for the models considered. Shen-EoS models are 
shown with a circle, APR-EoS models with triangles, and the 
LS-EoS model with a star. The two corotating models are 
indicated by ' x ', while the counterrotating models are marked 
with '-I-', The horizontal spread of the model group with the 
Shen-EoS is mainly caused by different total system masses. 
Lines connect cases with the same total mass and spin setup. 



little influence (Fig. [2l middle panels). Indeed, all mod- 
els using the Shen-EoS lead to values around 2.1-2.4 kHz 
for /peak J where the variation among the Shen- models of 
about 0.3 kHz is mostly due to the total system mass. On 
the other hand, the models using the APR-EoS and LS- 
EoS with their more compact NSs do not only yield larger 
values for /max [1, 0] but also much larger ones (around 
3.6 kHz) for /peak- The latter quantity mainly depends 
on the behavior of the EoS in the density regime between 
0.5po and 2po (po ~ 2.7 x lO^'^gcm"'^ being the nuclear 
saturation density), where the bulk of the remnant mass 
is located. While peak temperatures of several lOMeV 
are present in the HMNS at such densities (see 0]), T ^ 
contributions to the gas pressure affect the basic proper- 
ties of the post-merger oscillations only moderately. A 
comparison of models C1315, C1216, and A151414 with 
S1315, S1216, and A21414, respectively, shows that /peak 
decreases by at most ^0.2 kHz when thermal pressure is 
included (Table [l|. This is caused by the less compact 
structure of the HMNS in these cases. Considering the 
radiated energies before and after merging, Ai?i„ and 



Ai?pm, respectively, we find a characteristic variation of 
the GW signal with the NS spins. As shown in Fig. [3l 
the ratio of AiSpm to Ai?i„ is highest for counterrotating 
cases and lowest for corotating NSs. This is so because 
corotation leads to a stronger inspiral signal due to a pos- 
itive contribution from the NS spins, while damping the 
amplitude of the post-mergerpart due to a smaller non- 
axisymmetry of the remnant [7| . Counterrotation has the 
opposite effect (cf. Fig. [T]). The degeneracy of /max visi- 
ble in Fig. [3] for cases with APR and LS EoS can be lifted 
when the ratio Ai?pm/Ai?in is taken into account. 

A NS merger therefore produces a GW signal whose 
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FIG. 4: GW spectra 2f\h+{f)\ for two typical models with 
different EoSs, assuming a source distance of 20Mpc. Also 
shown is the strain noise amplitude /irms of advanced LIGO 
and of the resonant spheres detector DUAL. 



SNRs may be too low by up to a factor of 1.7. More- 
over, the post-merger signals are likely to be emitted for 
longer times than the considered window of 5 ms. Fitting 
an exponential decay to the GW amplitudes, we find a 
decay time of about 5 ms. From this, a further increase 
of the SNR by ~15% is estimated. Taking these cor- 
rections into account, we obtain for typical models like 
S1414 and APR21414 a SNR of -5 in LIGO and -6.5 
in DUAL. Assuming that a minimal value of about 3 is 
needed for detection in case the preceding inspiral chirp 
has been measured, such GW signals may be identified 
up to ~35Mpc (LIGO) and -45Mpc (DUAL). Accord- 
ing to 0, [3] , these distances correspond to event rates 
of 0.04-0.5 (LIGO) and 0.08-1.1 (DUAL) per year. 

We thank the DUAL research team for advice. This 
work was supported by the DFG through grants SFB- 
TR 7 and SFB-375. 



location in the space of the parameters /max, /peak, and 
Ai?pm/A£'in depends distinctively on the properties of 
the nuclear EoS. The latter determines the compactness 
of the merging stars and of the HMNS and thus the GW 
frequencies and energies emitted during the final plunge 
and post-merger oscillations. Characterizing a GW mea- 
surement in terms of the three parameters therefore pro- 
vides direct information about the NS EoS, in particular 
if the system mass is known from the inspiral chirp signal. 
This is a promising alternative to constraining the EoS 
by NS mass and very difficult radius determinations (e.g., 
[l9|). More work, however, is needed to understand the 
post-merger oscillations in terms of involved eigenmodes 
of the HMNS ([13 and refs. therein), and the GW pa- 
rameters for a large variety of EoSs should be computed 
by accurate GR merger simulations. 

To assess the detectability of the post-merger GW 
emission we follow Ref. [2l| and consider the ad- 
vanced LIGO interferometer and the omnidirectional 
DUAL detector, which consists of two nested reso- 
nant spheres [23] ■ The signal-to- noise ratio (SNR) for 
a given GW signal h{t) can be written as SNR^ = 
/^dln/|2//i(/)|V/irms(/), where h{f) is the fourier 
transform of h{t) and /irms(/) is the strain noise of the 
detector. In Fig. (H we take h{t) = h^{t) and compare 
the spectrum 2/|/i+| for two representative models with 
the strain noise of the LIGO and DUAL instruments. 
For an interferometer, h depends on the source orien- 
tation and direction relative to the interferometer arms 
and in the optimally aligned case is equal to /i+ or fi X ■ 
For DUAL, \h{f)\^ = \h{f)+? + \~Hf)x\^, because both 
polarisations can be measured simultaneously. In Ta- 
ble H] the SNRs are listed for all of our models. Note 
that only the post-merger waveforms are considered and 
a distance of 20 Mpc is assumed. Since we have underes- 
timated the wave amplitudes by ^40% (see above), our 
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